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The lepton flavour violating signals e + e~ — > £ + e~ and e~e~ — > £~e~ (£ — n,r) are studied 
in the context of low energy R-parity conserving supersymmetry at center of mass energies of 
interest for the next generation of linear colliders. Loop level amplitudes receive contributions from 
electroweak penguin and box diagrams involving sleptons and gauginos. Lepton flavour violation 
is due to off diagonal elements in SU (2) l doublet slepton mass matrix. These masses are treated 
as model independent free phenomenological parameters in order to discover regions in parameter 
space where the signal cross section may be observable. The results are compared with (a) the 
experimental bounds from the non-observation of rare radiative lepton decays /i, r — ■» ej and (b) the 
general mSUGRA theoretical scenario with seesaw mechanism where off diagonal slepton matrix 
entries are generated by renormalization group evolution of neutrino Yukawa couplings induced by 
the presence of new energy scales set by the heavy SU(2)l singlet neutrino masses. It is found that 
in e~e~ collisions the (er) signal can be observable with a total integrated luminosity of 100 fb _1 
and the that the background can be easily suppressed. In e + e~ collisions the cross section is smaller 
and higher luminosities are needed. The experimental bound on the decay n — > e-y prevents the 
(en) signal from being observable. 

PACS numbers: 11.30.Fs, 11.30.Pb, 12.gO.Jv, 14.80.Ly 



I. INTRODUCTION 

In the advent of growing evidence for neutrino oscil- 
lations and hence flavour mixing in the lepton sector of 
the Standard Model (SM), the topic of lepton flavour vi- 
olation (LFV) has received considerable attention. Non- 
vanishing neutrino masses in principle induce LFV pro- 
cesses such as £ — > I'j. If neutrinos have masses in the 
cV or sub-eV range, the neutrino generated branching ra- 
tio to the latter process is of order 0(1O~ 40 ) and there- 
fore unobservably small. On the other hand, in super- 
symmetric (SUSY) extensions of the SM the soft SUSY 
breaking potential V so ft contains, in general, nondiago- 
nal entries in generation space and therefore additional 
potential sources for LFV. Even in minimal supergravity 
scenarios characterized by universal soft mass term for 
scalar slepton and squark fields, renormalization induces 
potentially sizable weak scale flavour mixing [1] in V so ft- 

Many experimental efforts have been devoted to search 
for LFV and lepton number violating reactions, both 
in rare decays and in high energy accelerators. The 
strongest bounds on LFV come from the non-observation 
of radiative lepton decays [2-4] : 

Br(n -> 67) < 1.2 x lCT 11 , 

Br[T -> erf) < 2.7 x 1CT 6 , 

Br{r -> ni) < 1.1 x 1(T 6 . (1) 

The four LEP experiments searched for Z — > 
f = e,/i,T, i ^ j at the Z peak providing the 
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following upper bounds on branching ratios [5]: Br(Z — > 
e/x) < 1.7 x 1CT 6 , Br(Z -> er) < 9.8 x 1CT 6 , Br{Z -> 
t/i) < 1.2 x 10~ 5 . The high luminosity GigaZ option 
of the Tesla project [6] is expected to probe the above 
branching ratios down to ~ 0(1O -8 , 10 -9 ). A recent 
study of LFV induced by R-parity conserving SUSY at 
the Z peak is given in Ref . [7] . 

However it is interesting to know if such signals can 
be observed at higher energies. The OPAL collaboration 
searched for LFV reactions up to the highest center of 
mass energy reached by LEPII, y/s = 209 GeV [8]. One 
e + e~ — > efi event was found at ^fs — 189 GeV matching 
all tagging conditions, but it was interpreted as due to 
initial state radiation [8] . This negative result implies the 
following upper limits (at 95% confidence level) on the 
cross sections of LFV processes (for 200 GeV < *Js < 
209 GeV): 

cr(e+e- -> en) < 22 fb, 
a(e+e- -► er) < 78 fb, 
(j(e+e- -> nr) < 64 fb. (2) 

For limits corresponding to lower energies see Ref. [8] . 

In the following this approach will be pursued further 
and a detailed study of the reactions 

e+e- -► £+e~ , 

e~e~ -► t~e~ (£ = (J.,t) (3) 

will be presented in the context of SUSY extension of 
the SM with conserved R-parity for center of mass ener- 
gies of interest for the next linear collider projects (LC). 
The processes in Eq. (3) have the advantage of providing 
a clean final state being easy to identify experimentally 
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(two back to back different flavor leptons), though one 
has to pay the price of dealing with cross sections of 
order 0(a 4 ). Previous studies on SUSY induced LFV 
at a LC (see e.g. [9]) were limited to tree level pro- 
cesses for SUSY partner production decaying into final 
states characterized by very complicated topologies such 
as tilj + 4 jets + $ involving jets and missing energy. 
A detailed study of cuts and background is necessary to 
isolate the signal. 

The relevant Feynmann diagrams describing the pro- 
cesses in Eq. (3) are shown in Figs, f, 2, 3. They are 
the high energy analogue of the box and penguin dia- 
grams that mediate LFV rare decays as e.g. \i — > e + 7 
or fi — ► 3e. Due to the experimental limits on the cross 
sections and the loop nature of the process event rates 
are expected - even in more optimistic cases - to be rel- 
atively small. However, when the energy dependence of 
four-point and three-point functions is taken into account 
the amplitudes can show a resonance behavior as the en- 
ergy approaches thresholds for particle production. This 
is a consequence of the discontinuity of the derivative 
of the real part of a loop amplitude where it develops 
an imaginary part (Cutkosky rule). The cross section in 
this point may increase by orders of magnitude. We have 
shown in a recent paper [10] on LFV induced by heavy 
Majorana neutrinos that the enhancement may be quite 
dramatic in some regions of the parameter space. 

The plan of the paper is the following. Sec. II discusses 
LFV in R-parity conserving SUSY and gives an outline 
of the calculation. Sec. Ill contains numerical results 
for the signal cross section and a discussion of possible 
backgrounds. Sec. IV is devoted to a comparison with 
bounds from rare LFV lepton decays. Sec. V contains 
the conclusions. Appendices A and B give details of the 
lagrangians and numerical tools used in the calculation. 
Finally, in Appendix C hclicity amplitudes for e + e~ and 
e~e~ collisions are given. 

II. SUSY ORIGIN OF LEPTON FLAVOUR 
VIOLATION 

One of the most important challenges in contemporary 
particle physics is to understand the origin of neutrino 
masses. Quite generally this requires new fields to be 
added to those of the SM and/or those of its minimal 
SUSY version (MSSM). In the seesaw framework - the 
simplest scenario for the explanation of neutrino masses 
- and its SUSY extension, the superpotential contains 
three SU(2)l singlet neutrino supcrficlds Ni with the 
following couplings [11-15]: 

W = (Y v )ijE a H% Ni Lj + ^(Mn^NiNi. (4) 

Here H 2 is a Higgs doublet superfield, L t are the SU(2)l 
doublet lepton superfields, Y v is a Yukawa coupling ma- 
trix and Mr is the SU (2)l singlet neutrino mass matrix. 
As is usually done the basis has been chosen such that 



Mr is diagonal. The effective low energy neutrino mass 
matrix is given by 

M v = m T u M- y m u , (5) 

where mj = V2(Y„)ij /\/2 is the Dirac neutrino mass 
matrix and w 2 = (H®). 

Standard mSUGRA models contain an universal GUT 
scale (i.e. at the energy scale where the coupling con- 
stants unify) scalar field mass term mo. At low ener- 
gies the renormalization group equations (RGE) produce 
diagonal slcpton mass matrices. With the additional 
Yukawa couplings in Eq. (4) and a new mass scale (Mr) 
the RGE evolution of the parameters is modified: as- 
suming that Mr is the mass scale of heavy right-handed 
neutrinos, the RGE from GUT scale to Mr induce off- 
diagonal matrix elements in (m|)y. In the one loop ap- 
proximation the off-diagonal elements are [14]: 

«k * - g^(3 + al)ml(Y}Y v )a In (^g 1 ) ■ (6) 

Here an is a dimensionless parameter appearing in the 
matrix of trilinear mass terms At = Y^aoTOo contained in 

Vsoft ' 

The slepton mass eigenstates are obtained diagonaliz- 
ing the slepton mass matrices. The corresponding mix- 
ing matrices induce LFV couplings in the lcpton-slepton- 
gaugino vertices 'njihjX- The same effect on the 
mass matrix of SU(2)l singlet charged sleptons (m^)jj 
is smaller [14, 15]. 

The magnitude of LFV effects will depend on the RGE 
induced non diagonal entries and ultimately on the neu- 
trino Yukawa couplings (Y v )ij. These in turn depend on 
the fundamental theory in which this mechanism is em- 
bedded (for example SU{5) or 50(10) SUSY GUT [15- 
17]) and on the particular choice of texture for the neu- 
trino mass matrix [14, 18, 19]. The rate of LFV tran- 
sitions like £i — > £j, i ^ j, I = e,/j,,T induced by the 
lcpton-slepton-gaugino vertex is determined by the mix- 
ing matrix Uhij that, as stated above, is model depen- 
dent. In a model independent way, however, one can 
take the lepton, slepton, gaugino vertex flavour conserv- 
ing with the slepton in gauge eigenstates, so that LFV is 
given by mass insertion of non diagonal slepton propaga- 
tors [1, 7, 12]. 

In a similar spirit, the phenomenological study pre- 
sented in this paper will be quite model independent and 
in order to keep the discussion simple the mixing of only 
two generations is considered, so that the slepton and 
sneutrino mass matrix is: 

„ 2 ( m 2 Am 2 \ 

m L = y Am 2 ^ y (7) 

with eigenvalues: rh± = fn 2 ± Am 2 and maximal mixing 
matrix 
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Under these assumptions the LFV propagator in momen- 
tum space for a scalar line is 



(a) 



1 



1 



2 \^p 2 — mi_ p 2 — fh 2 _ 
Am 2 

(p 2 — m\ ) (p 2 — fh 2 _ ) ' 



(9) 



while a lepton flavour conserving (LFC) scalar line is de- 
scribed by 



2 \p z — m+ p z — m_ J 



(10) 



Therefore the essential parameter that controls the LFV 
signal is 



5ll = 



(11) 



Before presenting detailed calculations a qualititave order 
of magnitude estimate of the cross section can be given 
using dimensional arguments. Consider for simplicity a 
box diagram. Neglecting the external momenta in the 
loop propagators and indicating with rag a, typical SUSY 
mass, one has for the amplitude in the case of a scalar 
four point function: 



M 



9 



(4tt) 5 



(12) 



The constant comes from couplings and loop integration, 
the factor s from the spinorial part, the mass squared fac- 
tor from the numerator of the two gaugino propagators 
and the last factor from the loop integral. The corre- 
sponding total cross section (assuming polarized initial 
particles) is: 



1 

167T 



sin U\y 



S 2 — 

°LL 4 ' 

m% 



(13) 



Taking m s = 100 GeV, S LL = 0.1 and y/i = 200 GeV 
one has a ~ 1.3 X 10~ 2 fb while with y/s = 500 GeV 
a ~ 8 fb. With an annual integrated luminosity of order 
Lo=100 fb -1 one may expect an observable signal. 

However this estimate is clearly too crude: it gives a 
linear increase with s while one expects at high energies, 
y/s >> ms, a cross section which scales as a -1 . To get 
a realistic result it is necessary to compute exactly the 
energy dependence of the loop integrals and the interfer- 
ence among all contributing graphs. 



III. NUMERICAL RESULTS 

In the reactions considered here there are only leptons 
in the initial and final state. At the energies of a LC lep- 
ton masses can be safely neglected and thus all the cal- 
culations are done assuming massless external fermions. 
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FIG. 1: Box diagrams for e + e collisions. The full black dot 
in a scalar line denotes the lepton flavour violating propagator 
(Eq. 9). 



The signal is suppressed if neutralinos and charginos x '^ 
are Higgsino-like, since their coupling is proportional to 
the lepton masses. For the same reason left-right mix- 
ing in the slepton matrix is neglected. Therefore it is as- 
sumed that the two lightest neutralinos are pure Bino and 
pure Wino with masses M\ and Mi respectively, while 
charginos are pure charged Winos with mass M2, M\ and 
Mi being the gaugino masses in the soft breaking poten- 
tial. The relevant parts of the interaction lagrangian are 
listed in Appendix A. 

Due to the chiral nature of the couplings it is con- 
venient to calculate the amplitudes using the helicity 
base for spinors: the amplitudes are written in terms 
of spinor products and a numerical code can be easly 
implemented to compute both real and imaginary parts. 
Interference terms are also accounted for by summing the 
various contributions before taking the absolute modulus 
squared of the amplitude. In the helicity basis and in the 
limit of massless fermions there are only two independent 
spinors, u+(p) = Ur(p) and = Ul(jp) with only 

two non-zero spinor products: un(pa)uL(pb) = S(p a ,Pb), 
u L {pa)uB{Pb) = T(p ai p b ) given by compact expressions, 
see Appendix B 1. The loop integrals are decomposed in 
form factors and calculated numerically using the pack- 
age LoopTools [20] . The decomposition of loop inte- 
grals is obtained for massless external particles and with 
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FIG. 2: Penguin and external legs diagrams for e + e~ colli- 
sions. The full black dot in a scalar line denotes the lepton 
flavour violating propagator. In the diagrams where it is not 
marked it can occour in both lines. Diagrams like (d) and (e) 
with the gauge boson in the s channel are also present but 
not shown. 



the loop momenta assigned as described in Fig. 12 of 
Appendix B 2. The exact dependence from the masses of 
the particles exchanged in the loop is also given in Ap- 
pendix B 2. Assigning the momenta in a different way 
corresponds to a shift of the integration variables and 
produces different combinations of the loop form factors 
appearing in the amplitudes. The numerical values re- 
main unchanged. 

Besides computational advantages the helicity method 
clarifies the physics of the processes. The momenta of the 
external particles are specified as in Eq. (B3) (Appendix 
B) and Fig. 12 (Appendix C), and the following reactions 
are considered: 



Oi,Ai)e (p 2 ,A 2 ) 
"(Pij Ai)e"(p 2 , A 2 ) 



(■ (P3, A 3 )e + (p 4 , A 4 ), 
rCp3,A 3 )e-(p4,A4). (14) 



Here Ai denotes the helicity of particle i. The correspond- 
ing helicity amplitudes Mj expressed in terms of spinor 
products and LoopTools form factors are obtained af- 
ter tedious but straightforward algebra. They can be 
found in Appendix C. 

The integrated cross sections corresponding to each in- 



FIG. 3: Feynman diagrams for e~e~ collisions. The full black 
dot in a scalar line denotes again the lepton flavour violating 
propagator. Exchange diagrams are not shown. 



dividual amplitude Mj is: 
1 



327TS 



d(cos6) \Mj 



(15) 



The total unpolarized cross-section (averaged over ini- 
tial spins) is a = (l/4)X)j cr j- The dependence on 
the scattering angle is encoded in the Mandelstam vari- 
ables u and t. Numerical results are obtained using 
the mSUGRA relation M\ ~ O.5M2 for gaugino masses 
while Am 2 and the slepton masses are taken to be free 
phenomenological parameters. The parameter space is 
scanned in order to identify the regions which may de- 
liver an interesting signal. The discussion of whether such 
regions are compatible with present experimental bounds 
is postponed to the next section. 



The contributing amplitudes are (Appendix CI): 



M A 
Mb 
M c 



M{e 
M{e 
M(e 
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r L e + L ). 



(16) 
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FIG. 4: Differential cross section as a function of the scat- 
tering angle for e + e~ collisions. The following values of the 



parameters are used: Mi = 80, 
GeV and Am 2 = 6000 GeV 2 . 



M 2 



160, m~ t — mo 



100 



For each helicity amplitude the corresponding differential 
polarized cross section is shown in Fig. 4. The different 
behavior is easily understood in terms of helicity conser- 
vation at high energy. Amplitude Ma is peaked in the 
forward direction since it has a P-wave initial state with 
J z = +1. Angular momentum conservation requires the 
right-handed positron to be emitted in the positive direc- 
tion of the collision axis while the left-handed negative 
charged lepton must have its momentum in the opposite 
direction. Amplitude Mb is peaked in the backward di- 
rection as it is a P-wave scattering with J z = — 1. The 
right-handed positron must be emitted backward while 
the negative charged lepton is in the forward direction. 
Amplitude Mc has no virtual vector boson exchanged 
and is an S-wave (J z = 0) scattering. One expects there- 
fore an almost flat, isotropic distribution. 

The dominating contribution to the integrated unpo- 
larized cross section comes from amplitude Ma, which is 
an order of magnitude larger than M b and two orders of 
magnitude larger than Mc ^ most of the phase space. 
Only for large scattering angles (backward direction) the 
amplitude Mb dominates and Ma is the smallest one. 
In Fig. 4 the dotted line corresponds to the unpolarized 
differential cross section (i.e. the inchocrent sum of the 
contributions of Ma,b.c averaged over the initial spins). 
It is worth to remark that in such circumstances the pos- 
sibility of having polarized electron and positron beams 
would maximize the chances to observe these signals. 
Considering the unpolarized cross section corresponds es- 




FIG. 5: Total cross section (fb) as a function of ^fs for e + e 
collisions for the three helicity amplitude. The parameters 
chosen are Mi = 80, M 2 = 160, m~ t = mo = 100 GeV. The 
solid lines correspond to Am 2 increasing from 100 GeV 2 to 
900 GeV 2 in steps of 100. The dashed lines correspond to 
Am 2 increasing from 1000 to 8000 GeV 2 in steps of 1000. 



sentially to calculating a unpo i ~ (1/4) cr(e^e£ — > ^£e^). 

Fig. 5 shows the cross section integrated over the scat- 
tering angle for the three helicity amplitude as a function 
of the center of mass energy y/s and for increasing values 
of the LFV parameter Am 2 . The presence of spikes is 
due to the onset of the absorptive part of the diagrams 
corresponding to thresholds of real particle pair produc- 
tion. For the values of masses used in Fig. 5 one expects 
thresholds effects at <~ 200 GeV for slepton pair produc- 
tion and ~ 320 GeV for gaugino pair production. This 
is evident for in a a (upper- left panel) and <jb (upper- 
right panel). The shape is determined in the first case 
by the destructive interference among the two types of 
box graphs (with scalars and fermions on threshold) and 
by the value of Am 2 inducing two distinct thresholds at 
to 2 ± Am 2 , ob is determined only by penguin diagrams 
that give smaller contribution relative to the boxes, ac 
receives contributions only from box diagrams: at thresh- 
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FIG. 6: Total cross section for the amplitude A as a function 
of the dimensionless parameter 8ll (see Eq. (11)) and for 
yfs = 2mi. The values of the other parameters are given in 
the legends. 



old for sleptons production its value varies by orders of 
magnitude differently from the two other cases. This can 
be easily understood considering the threshold behav- 
ior of the cross section for sleptons pair production [21]: 
defining (i the selectron velocity, the intermediate state 
of the amplitudes Ma and Mb correspond to the reac- 
tions e^e^ — > e£ej and e R e^ — > e^ej that near thresh- 
old behave like /3 3 while amplitude Mc corresponds to 
the reaction e L e\ — > e L e~^ that at threshold behaves like 
P. 

The cross section is peaked around y/s = 2rfiL- In 
Fig. 6 <ta = a{e R eJ j — > -^ e j?) is shown as a function of 
5ll for V* = 2tol and for different values of slepton and 
gaugino masses. Given an annual integrated luminosity 
L =100 a cross section of 10~ 2 fb produces one sig- 
nal event per year. Such an event rate is reached only for 
Mi not larger than - 200 GeV and 6 el — 0(1). This hy- 
pothesis will be discussed in the next section. Moreover 
angular cuts in the forward direction are needed to sup- 
press possible SM backgrounds and - since the largest 
values of the cross section correspond to small angles - 
the signal will be affected by such a cut. Therefore the 
observation of LFV in e + e~ collisons appears to be diffi- 
cult unless L is considerably larger than 10 2 fb _1 /yr. 



B. e e collisions 

The contributing amplitudes are (Appendix C 2) 

M E i = M(elel^eiel), 
M E 2 = M(e L e R -> l L e R ), 
M E3 = M{e R e L ^ l L e R ). (17) 

The corresponding differential cross sections are plot- 
ted in Fig. 7. Mei has J z = 0, is flat and forward- 
backward symmetric because of the antisymmtrizzation. 
M E2 and M £3 describe P-wave scattering with J z = +1 
and J z = — 1 respectively: in order to conserve angular 
momentum M E2 must be peaked in the forward direction 
while M E3 favours backward scattering. Both M E2 and 
Me3 are orders of magnitude smaller than Mei- The 
signal cross section is to a very good approximation given 
by the amplitude Mei- Since it is almost flat the angu- 
lar integration will give a factor almost exactly equal to 
two. This again shows the importance of the option of 
having polarized beams. If both colliding electrons are 
left-handed one singles out the dominant helicity ampli- 
tude and a factor four is gained in the cross section rel- 
ative to the unpolarized case. This may be important 
in view of the relatively small signal cross section one is 
dealing with. In this case, due to the smaller number 
of diagrams, the analysis of the total cross section as a 
function of y/s is easier (see Fig. 8): the box diagrams 
dominate at y/s — 2m e where a changes of orders of mag- 
nitude giving a sharp peak that is smeared only by large 
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FIG. 8: Total cross section of the dominant amplitude El 
as a function of yfs for e~e~ collisions. The choice of the 
parameters is the same as in Fig. 4. 



values of Am 2 , while penguin diagrams give a substan- 
tial contribution only at higher energies. The reason is 
the same as for the ac behaviour in the e + e _ case: the 
intermediate state e^e^ — > e^e^ behaves like (3, while 
the other two like (3 3 . Here the highest absolute value 
is due to the couplings and the constructive interference 
of boxes where both Binos and Winos can be exchanged. 
The dependence of gei on Sll is shown in Fig. 9. With 
SUSY masses not much larger than ~ 200 GeV the signal 
is of order 0(1Q- 2 ) fb for 5 LL > ©(lO" 1 ). Relative to 
the e + e~ case there are two important features: (i) the 
cross section is practically angle independent so that it 
is insensitive to angular (or tranverse momentum) cuts; 
(ii) the SM background - though not completely absent 
- can be easily controlled as will be shown in the next 
subsection. 



1. Background 

The signal has the unique characteristic of a back to 
back high energy lepton pair. Sources of background were 
qualitatively discussed in Ref. [22]. 

Initial and final state radiation can be a source of back- 
ground. An example is the OPAL event [8], although lep- 
ton pairs can hardly have the same kinematical feature 
of the signal. Other sources present multiparticle final 
states (at least six particles) and missing energy due to 
the presence of neutrino pairs. 

The first type is given by reactions like e~e~ — > 



o P1 (fb) 




FIG. 9: Total cross section for the amplitude El in function 
of the dimensionless parameter Sll, see Eq. (11). The values 
of the other parameters are given in the legends. Each plotted 
line is calculated assuming y/s = 2rh l ■ 



e~e~b*b* that proceedes trough virtual photon fusion. 
The subsequent chain of weak deacys produces a final 
state with missing momentum, hadronic jets and oppo- 
site or same sign leptons, that however can be again sep- 
arated using the clear kinematical topology of the signal. 
A second type: 



e e — > u e v e W W 

W~*W~* -» trvitr 



(18) 



with four neutrinos and a like sign-dilepton pair that can 
be of the same or different flavour. This appears to be the 
most dangerous background, as it produces two leptons 
and missing energy, and therefore it is analyzed in more 
detail. Moreover, to the best of our knowledge, it has 
not been previously considered in the literature. Fig. 10 
shows the total cross section e~ e~ — > v e v e W~W~ cal- 
culated with the CompHEP package [23], that allows 
to compute numerically the 17 Feynman diagrams con- 
tributing at tree level. Above the threshold for W~W~ 
gauge boson production the cross section rises rapidly by 
orders of magnitude, becoming almost constant at high 
energies. In the region y/s ~ 250 — 400 GeV it increases 
from 10~ 2 fbtolfb. In order to get an estimate of the 
cross section for the six particle final state process, the 
cross section a(e~e~ — > W~W~ vv) has to be multiplied 
by the branching ratio of the leptonic decays of the two 
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FIG. 10: Total cross section and distributions for e~e~ — > 
W~W~h>v Upper-left figure: total cross section as a function 
of y/s. Upper-right: angular distribution for a W~ where 9 is 
the angle among the collision axis and the boson momentum. 
Bottom-left: distribution of the transverse momentum of W~ . 
Bottom-right: energy distribution of the two neutrinos. All 
distributions are calculated with y/s = 300 GeV. 



gauge bosons, ~ 10%, so that cr B ackground — 10 4 - 10 2 
fb, and it is at the level of the signal. However the kine- 
matical configuration of the final state leptons is com- 
pletely different. Fig. 10 (upper-right) shows the angu- 
lar distribution of the gauge bosons which is peaked in 
the forward and backward directions so that the leptons 
produced in the W gauge boson decay are emitted pref- 
erentially along the collsion axis. In addition their tans- 
verse momenta will be softer compared to that of the 
signal: Fig. 10 (bottom-left panel) shows that the trans- 
verse momenta distribution of the gauge bosons is peaked 
at p£ = (yfs/2 ~ M w )/2 ~ 35 GeV for yfs = 300 GeV. 
Consequentely the lepton distributions will be peaked at 
Pt/2 — 17.5 GeV. The missing energy due to the unde- 
tected neutrinos (Fig. 10, bottom-right panel) can be as 
large as ~ \/s — 2M\y- This distribution should be con- 
voluted with that of the neutrinos produced in the gauge 
boson decay. Therefore it can be safely concluded that it 
will be possible to control this background because, with 
reasonable cuts on the transverse momenta and missing 
energy, it will be drastically reduced, while - as men- 
tioned above - these cuts will not affect significantly the 
signal. 



IV. COMPARISON WITH RARE LEPTON 
RADIATIVE DECAYS 

The main result of the calculations presented in the 
previous sections is that, as it can be inferred from 
Fig. 9, the phenomenological points of the SUSY param- 
eter space corresponding to gaugino masses (Mi, M 2 ) = 
(80, 160) GeV or (100, 200) GeV and to slepton masses 
m L = 100 - 200 GeV and S LL > lO" 1 (which implies 
Am 2 > 10 3 GeV 2 ) can give in the e~ e~ mode a de- 
tectable LFV signal (e~e~ — > i~e~) although at the 
level of 0(1 - 25) events/yr with L = 100 fb" 1 . Higher 
sensitivity to the SUSY parameter space could be ob- 
tained with larger L - It is interesting to note that this 
light sparticles spectrum that is promising for collider 
discovery, is also preferred by the electroweak data fit. 
In Ref. [24] it is shown that light sneutrinos, charged 
left sleptons and light gauginos improve the agreement 
among the electroweak precision measurements and the 
lower bounds on the Higgs mass. 

On the other hand the experimental bounds on rare 
lepton decays set constraints on the LFV violating 
paramters Am 2 or Sll- the constraints in Eq. (1) de- 
fine an allowed (and an excluded) region in the plane 
(SlLtItt-l) which are computed using the formulas given 
in Ref. [13] (adapted to our model) for the LFV radiative 
lepton decays. These regions have to be compared with 
those satisfying the "discovery" condition 

L a(5 LL ,m L )>l. (19) 

Such a comparison is shown in Fig. 11 from which it 
emerges that: (i) For the e~e~ — > £~e~ process there 
is an observable signal in the upper left corner of the 
(Sll, mi) plane. The extension of this region depends on 
Lq. (ii) The bound from r — > ej does not constrain the 
region of the (Sll, "II) plane compatible with an observ- 
able LFV signal and therefore the reaction e~e~ — > r~e~ 
could produce a detectable signal whithin the highlighted 
regions of the parameter space (upper-left regions in the 
(Sll, mi.) plane), (in) As regards the constraints from 
the /j, — > e-y decay the allowed region in the (Sll^l) 
plane is shown by the circular dark dots (red with colour) : 
the process e~e~ — ► /i~e - is observable only in a small 
section of the parameter space since the allowed region 
from the \i — > decay almost does not overlap with the 
collider "discovery" region except for a very small frac- 
tion in the case of gaugino masses (Mi = 80 GeV and 
M2 = 160 GeV). The compatibility of values of Sll ~ 1 
is due to a cancellation among the diagrams that describe 
the i — > £'7 decay in particular points of the parameter 
space [7]. 

As regards the radiative mechanism that generates the 
off-diagonal elements in mSUGRA models (as discussed 
in Sec. II) one should check if this mechanism may gen- 
erate large values of Am 2 . The answer is yes, at least 
for some particular scenario of neutrino masses and mix- 
ing. It is well known that any 'bottom-up' approach that 
reconstructs the Y v from the see-saw mechanism and neu- 
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M J= 80 GeV, M 2 = 160 GeV 




-3 I* I I * I i i ■ i I *ii 

1U LOO 200 300 400 500 600 

(GeV) 

FIG. 11: Scatter plot in the plane {5t,L,m£) of: (a) the ex- 
perimental bounds from /j, — » e*y and r — » fi-y (allowed regions 
with circular dots); (b) regions where the signal can give at 
least one event with two different values of integrated lumi- 
nosity (squared dots), for two sets of gaugino masses. Each 
signal point is calculated at ^fs = 2rfiL- 



trino masses and mixings is ambiguous up to a complex, 
orthogonal matrix R [14]. Usually this matrix is taken 
to be real or identical to the unit matrix. However in 
Ref. [25] it is shown that in the case of a quasi-degenerate 
neutrino mass spectrum R being complex allows for val- 
ues of Am 2 being larger by 5-8 orders of magnitude rel- 
ative to the case of R being real or the unit matrix. In 
this case one has [25] 



V. SUMMERY AND CONCLUSIONS 

The search at lepton colliders of lepton flavour and 
lepton number violating signal is complementary to the 
search for rare leptons decays. The next generation of 
linear colliders will offer the opportunity to look for reac- 
tions like e ± e~ — > ^ ± e _ (£ = fi,r) at energies well above 
the Z peak resonance. Upper bounds on the cross sec- 
tions for these processes at the highest energies reached 
by LEP, 189 GeV < y/s < 209 GeV, where given by the 
OPAL collaboration, Eq. (2). 

In this paper the reactions e^e - — > £^e~ (£ — ^,r) 
induced by sleptons mixing in R-parity conserving su- 
persymmetry have been studied. The reactions proceed 
through loop diagrams (box and penguin type) involv- 
ing sleptons, neutralinos and charginos. The amplitudes 
have been evaluated in the helicity basis and the loop 
integrals are calculated numerically. The resulting cross 
sections exhibit the well known threshold enhancement 
for center of mass energies corresponding to pair produc- 
tion of supersymmetric particles. In particular due to the 
dominance of the (s, i)-channel box diagrams with slep- 
tons on the threshold in the intermediate state, the LFV 
cross section reaches its maximum value at the energy 
corresponding to the threshold for sleptons pair produc- 
tion both in e + e _ and e~e~ collisions. 

The e~e~ option with left-polarized beams stands bet- 
ter chances to provide a detectable signal. A comparison 
with present experimental bounds on radiative lepton de- 
cays shows that an observable (e~ e~ — > r~e~) signal is 
compatible with the non observation of the decay r — > e"f 
giving some tens of events with an integrated luminos- 
ity of 100 fb _1 . On the contrary the more restrictive 
constraints from the non-observation of fi — > make 
the search of e~e - — > fi~e~ unrealistic unless the inte- 
grated luminosity is very large. It has been shown that 
the Standard Model background is low and can be easily 
suppressed using that the signal final state consists of two 
back to back high energy leptons of different flavour with 
no missing energy. The observation e + e~ in collisions 
will be more difficult because of smaller cross sections. 
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\(YjY„) 13 \ 2 ~ x O(0.1 - 1.0). (20) 



Choosing for example M R = 2-10 14 GeV, M GUT = 2-10 16 
GeV, m v = 0.3 eV, v = 174 GeV in Eq. (20) and a = 
0, m = 150 GeV in Eq. (6) Ato 2 varies in the range 
2400 - 7800 GeV 2 , i.e. with 100 GeV < m L < 200 GeV, 
8ll is in the range (0.06 — 0.78). 



APPENDIX A: LAGRANGIAN AND COUPLINGS 

The interaction lagrangians in the gauge basis for su- 
perparticles in the notation of Ref. [26] : 
(a) Lepton- chargino-sneutrino: 

C = 0flP R W9 + h.c, (Al) 

with coupling = —g. 
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(b) Lepton-neutralino-slepton: 

C = (0f 3 lP R W 3 L + 0f£P R BL + 0pP L BR) 



+ h.c, 



(A2) 



with couplings given by Of 3 = Of = -^t w , Of = 

V2gtw and L = £l, R = £ R . 
(c) Lepton-lepton-vector boson: 

£= E V~£^{O^P L + O^P R )l (A3) 



where Of = -gs w tw, O|o = +^7(5 - s w)-< ° 
0* = e. 

(d) Slepton-slepton-vector boson: 



C = iO^V^L* d "L, 



with — p — rt-S^fi — "l — 2_ 

(e) Chargino-chargino-vector boson: 



with 



-e, Ojo = -gc w - 



(A4) 



(A5) 



The spinor products are determined by the components 
of these four momenta in the following way: 



S(Pa,Pb) 



ijPl+ipl) 
{pl+ipl) 




V ' 

pi- vl 



(B4) 



and T(p ai p b ) is easily deduced by ralations (Bl). Using 
Eq. (B3) and Eq. (B4) it is easy to see that the rela- 
tions (Bl) are satisfied. In the case of 2 — > 2 scattering, 
with the momenta given in Eq. (B3), the preceding ex- 
pressions simplifies to: 

S{p a ,Pb) = (p Z a -Pb)+i(Pa-Pb), 

T(p a ,p b ) = (P z b-p z a )-i(p x b -p x a), (B5) 

and product of spinor products are directly related to 
s,t,u. For example one has: 



S( Pl ,p 3 )T(p 4 ,p 2 ) = 

S{ Pl ,P4)T{p S ,P2) = 

S(p 1 ,p 2 )T(p i ,p 3 ) = se 



-(1 + cos6») = u, 
|(l-cos0)=t, 



(B6) 



APPENDIX B: NUMERICAL TOOLS 
1. Spinor products 

Here the basic formulas used in the computation of 
helicity amplitudes are given. More details and proofs are 
given in Ref. [27]. The spinor products satisfy exchange 
relations: 

S(p a ,Pb) = S(p b ,Pa), T(p a ,p b ) = -T(p b ,p a ), 
S(p a ,p b )=T*(p b ,p a ), T(p a ,p b ) = S*(p a ,p b ), (Bl) 

\S{pb,Pa)\ 2 = 2pa-Pb, \T( Pa , Pb )\ 2 = 2p a - Pb . 

The necessary relations to write the amplitudes in terms 
of spinor products are the Chisholm identities: 

[uxipah^Uxiph^Jf, = 
2 [u x (p b )u\(p a ) + U-x{p a )U-xiPb)] , 

ft = u R (p)u R {p) + u L (p)u L (p). (B2) 

where A = L,R indicates the helicity of the spinor. The 
external momenta are parametrized in terms of the Man- 
dclstam variable s and the scattering angle in the center 
of mass frame: 



(1,0,0,1), 

(1,0,0,-1), 

(1, — sin 6*, 0, — cos#), 

(1, sin 0,0, cos 6). (B3) 





V~s 


Pi = 


2 




V~s 


P2 = 


2 




^ 


P3 = 


2 




yfs 


Pi = 


2 



2. Tensor integral decomposition 

The loop integrals are evaluated numerically with the 
package LoopTools [20] . Here we report the definitions 
and the decomposition for two, three and four point ten- 
sor functions. 



B 



d 4 q 



'' ~ 1 in 2 N 1 N 2 ) 

d 4 q q^qaqp 



C^-afj — - J 



D 



in 2 N!N2N 3 ' 
d 4 q q^qaqp 



(B7) 



J in 2 Ni^^Ni' 
are expressed as: 

Bft = ki^Bi 
2 

i=l 

2 

Cfiu = g^v Coo + k ^ k 3v C i3 

3 

Dfi = y ' ki^Di, 

i=l 

3 

= g^Doo + E ki^kj^ij, (B8) 

where the kiS are sums of external momenta appearing 
in the loops propagators as reported in Fig. 12: 

qi = q + ki = q + pi, 
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i>2- 



(a) 

<h 



(b) 



Pi- 



}> 2 . 



P\ 



Pi- 



9i 




93 




9 

(c) 






92»,i 




91^^93 



-P4 



(d) 



(e) 

a,' 



P3 P2 



P4 Pi 



-P3 P'l- 




(f) 9 



Pl- 



-P4 Pi. 




Pi 



-Pi 



FIG. 12: Definition of virtual momenta for kinematics and 
tensor integral decomposition. 



92 


= q 


-l-fca = 


9 + 


PH 


"P2, 


<Hu,t 


= q 


f Pl- 


P3, 






93 


= q 


f *3 = 


9 + 


P4, 




93 „ 


= q 


fP3, 








9i 


= q 


f fci = 


9 + 


P4, 




i 

<h 


= q 


f 4 = 


9 + 


P3 ^ 


"P4, 




= q 












= q 


fP3 - 


P2, 






q' 


= q 


-P3 









(B9) 



and the mass and Mandelstam variables dependence for 
a generic two, three, four point function and for the vari- 
ous topologies of graphs corresponding to the kinematical 
channels are: 

D a = D(0,0,0,0,S,t,m^, 771^,771^,771^), 

D b = D(0,0,0,0,s,u,m 2 ,m 2 qi ,m 2 q2 ,m 2 3 ), 
D c = D(0,0,0,0,u,t,m 2 q ,m 2 qi ,m 2 q2 ,m 2 q3 ), 

C d = C(0,0,S,T77y, 771^,771^), 

C e = C(0,0,t,m 2 q ,,m 2 qli ,m 2 q , 2 ), 



B f = B(m ,m q ,). 



(BIO) 



APPENDIX C: HELICITY AMPLITUDES 
1. e+e - collisions 

The amplitudes are given assuming that the negative 
charged final leptons has changed flavour. The other pos- 



sibility is taken into account simply multiplyng the total 
cross section by two. The non-zero helicity amplitudes 
are found to be: 



A: e+el 



L R 



For clarity, graphs are grouped according to the virtual 
particles present in the boxes that can be produced in 
e + e~ collisions: 
1. Virtual selectrons pair: 

There are four box diagrams with all the possible B and 
W assignament in the neutralino lines in Figs, la, lb: 



M D Ail = 



i,i'=B,W° 



(o l Y(o l YT(p 1 ,p 3 )s( Pi ,p 2 ) 



x {2D$(*,t) +T(pi,p 4 .)Sfa,p 1 )D$ (s,t) 
- mimi>iy$\u,t)}. (CI) 

The terms depending on (u, t) come from the crossed 
box diagrams due to the Majorana nature of neutralinos. 
Contribution from s and t channel penguins, Figs. 2a, 2b 
and the corresponding external legs corrections, Fig. 2d, 
give: 

M ai = 2T(pi,p 3 )5'(P4,P2) 



V,i 



(O i f(OY) 2 (B +B 1 )] 
Y J Dv{t)\{O l ) 2 O v "o v L 2C m (t) 



V,i 



(O i ) 2 (OY) 2 (B +B 1 )}} 



(C2) 



The photon and Z° propagators are given by Dy(s) = 

-i/(s - Ml + iM v T v ) (V = ~f,Z) for the s-channel, 

while no immaginary part is present in the denominator 

for t and u channels. 

2. Virtual sneutrinos pair: 

The box diagrams in Fig. lc reads: 



(O^) i T{p 1 ,p 3 )S< J)4 ,p 2 ){2D 00 (s,t) 



TfrupJSfapjDufrt)}, (C3) 
while the penguins of Figs. 2b, 2e: 

•^,4,2 = 2T (Pi:P3)S , (p 4 ,P2) 

x [d z (s) [{of) 2 OfO z L 2C m {s) 

(0t) 2 (0 z L ) 2 (B Q + B 1 ) 

D z (t) [(of) 2 OfO z L 2C 00 (t) 

- (Of) 2 (0f) 2 (i?o + i?i)]} (C4) 

The amplitudes present the same structure of those in 
case 1. 

3. Virtual chargino pair: 
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The box diagram in Fig. Id: 



M 



M% 3 = (of) i T(p 1 ,p 3 )S(p 2 ,p 4 ){2D oa ( S ,t) 

+ T( Pl ,p 2 )S( Pl ,p 2 )V(s,t)} . (C5) 

where k,l — s,t,u and 

V(k, I) = [D 12 (k, I) + D 22 (k, I) + D 23 (k, I) + D 2 (k, I)}, 
The penguin diagrams in Figs. 2c, 2e give: 



A 
B,l 



= 2T(p 1 ,p 4 )S(p 3 ,p 2 ) 

[ V.i 

- (0*) 2 (0}0 2 (5o + £i)] 
- ^^y(t)[(O 4 ) 2 O^2C 00 (t) 

{af{O v R f{B, + B 1 )}} 



(C8) 



M% 3 = 2T(pi,P3)S(p4,P2) 

x I E *m«) [(<# ) 2 <orc( S ) 

- (Of) 2 (0l) 2 (i?o + i?i/ 

- ]T D v {t) [(Of f0f O v L C{t) 

V=i,Z a 

- (OffiOlfiBo + B,)]} (C6) 
where k — s,t,u and 



C(fc) = {Co(fc) 



?71 ~ 

- [2C 00 (fc) + fc(C 2 (fc) + C 12 (k) + C 22 (k))}} . 
4- Virtual neutralino pair: 

There are four box diagrams with left-slepton and all 
possible combinations of Bino and neutral Wino in the 
loop of Figs, le, If with left sleptons exchanged: 



M D A,4 = ^(O t ) 2 (0 4 ') 2 T(p 1 ,p 3 )5b 2 ,P4) 

x +T(p 1 ,p 2 )S(p 1 ,p 2 )V i ' i '(s,t) 



m;m;>D\ 



?'(*>«)} 



(C7) 



Note that there is no penguin contribution to this 
channel in the B, W° basis. The amplitudes 3 and 
Ai 1 ^ 4 have a minus sign relative to the other ampli- 
tudes, because T(p 1 ,p 3 )S(p 2 ,p 4 ) = -T(p 1 ,p 3 )S(p 4 ,p 2 ), 
see Eq. (Bl). Its origin is due to the fact that once one 
fixes the order of the spinors, the two different topologies 
of box diagrams need an odd number of permutation of 
fcrmion fields to bring them in the same order. The same 
holds for the relative sign between s and t channel pen- 
guins. 



B: e\e R 



^L e R 



This differs from the previous by the exchange of initial 
state helicity: only the penguin diagrams contribute and 
the amplitudes are obtained selecting the O r ,Z P R oper- 
ator in the lepton-lepton-vector boson vertex. 



M B,2 = 2T (Pi^P4)S(p 3} p 2 ) 

x {d z (s) [(0ff0f0 z R 2C m ( s ) 

(Off{O z R )\B + B 1 ) 

- D z (t) [(0ff0f0 z R 2C w (t) 

- (Of) 2 (0|) 2 (5o + Si)]} 



M% z = 2T( Pl ,p 4 )S(p 3 ,p 2 ) 

x I E D v { s )\(0tf0f0 v R C( S ) 



V=i,Z° 



(O?y(O R y(B + B 1 )_ 
- Yl D v (t)\(0tf0t0 v R C(t) 



V=~i,Z° 



(C9) 



(otf(O R f(B + Bl )\} 



(CIO) 



C: e+el 



The box diagram in Figs, le, If with the right-handed 
selectron and the Binos in the neutralinos lines con- 
tributes: 



Mc,i = (Of) 2 (ol) 2 {S(p u p 2 )T(p 3 ,p 4 )2D oa ( S ,t) 

+ S{pi,p 2 )T{p 2 ,p 4 )T{p 3 , P i)S{pi, P2 )T>(s,t) 

- S(pi,p 2 )T(p 3 ,p 4 )[2D o(s,u) 

+ S(p l7 p 4 )T(p u p 4 )D 13 ( S ,u)}} (Cll) 

The box diagrams in Figs, lg, lh: 



M a Ct2 = (of) 2 (of) 2 s( Pl ,p 2 )T(p 3 , P4 ) 

x {m|A>(M) -2D Q0 (u,t) 

+ S( Pl , P4 )T( Pl , P4 )D 13 (u,t)} (C12) 



e e collisions 



The helicity amplitudes are: 



E1: e L e L -^^L e L 
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Four box diagrams of the kind given in Fig. 3a with left 
sleptons and the box in Fig. 3b with charginos: 



M D E1 = ^(0 4 ) 2 (O l ') 2 5(p 1 ,p 2 )T(p 4 ,p 3 ) 

x {mimi>Di{s,t) + 2D%(u,t) 

- S(p 4 ,pi)T(p 4 ,pi)D^ («,*)} 

+ {Off [2D c 00 (u,t) + S(p 4 , Pl )T(p 4 , Pl )D c 13 (u,t)} 

- M(pi^p2,u^t) (C13) 

Penguin diagrams in t and u channels, with left cou- 
plings with gauge bosons: 

M£ M - 2S( P 2,Pi)T(p 4 ,p 3 ) 

J2Dv(t) [{O l fO^O v L 2C m {t) 

V.i 



(OriOinBo + B,)} 
Y^Dviu) \(0*fO^O v L 2C m {u) 



V,i 



(O'fiOlfiBo + B,)]} 
M{p\^p2,u^t) (C14) 



M% 1<2 = 2S(p2,Pi)T( P4 ,p 3 ) 

x {D v (t) [{Of) 2 O™O z L 2C 00 (t) 

- (OriOinBo + B,)} 
- D z (u)[{of) 2 OfO z L 2C m (u) 

(O?f(O z L f(B + Bl )]} 

- M(pi^p2,u^t) (C15) 



M D E2 = (Of) 2 (0?) 2 T( Pl7 p 3 )S( P4 ,p2) 



x [2D 00 {s,t)+T(p 1 ,p 4 )S( P4 ,p 2 )D 31 {s,t) 

2 "./.. *VI (C17) 



■^£2,1 = 2T (Pl )P3) S(p 4 ,P2) 



^^(t){(O l ) 2 O^2C 00 (t) 
(O l ) 2 (0^) 2 (So + Si)} 



(C18) 



•^£2,2 = 2T (P1,P3)S(P4,P2) 

x J D z (t){(of) 2 Orol2C 00 (t) 

(of) 2 (Ol) 2 (B + Si)} (C19) 



■^£2,3 = 2T (PliP3)S(p4,P2) 

x ^ A,(t){(0f) 2 0^C(t) 

- (of) 2 (O^) 2 (B + Si)} (C20) 

This is obtained simply exchanging pi <-> j»2 and t <-» u 
in the previous amplitudes. 



- (Of) 2 (0|) 2 T( P2 ,p 3 )5(p4,Pl) 

x [2Doo(s,u) +r(p 2 ,P4)5'(P4,Pi)£ ) 3i(s,u) 

2 „.M ( C2 l) 



m~j£>o(t,u)] 



A 

£1,3 



M% 31 = 2T(p 2 , P3 )^(P4,Pi) 



= 2S( P 2,Pi)T{p 4 ,p 3 ) 

x I J] D v (t) [(Of fof O v L C(t) 

- (ofnoinBo+B,) 

- 53 JM«) [(Of) 2 O^OrC( U ) 

V=7,Z° 

- (of) 2 (OZ) 2 (So + Si)]} 

- M(pi^p2,u^t) (C16) 



All amplitudes are anty-symmetrized respect to initial 
state identical leptons. 



E2: e L e R 



Box diagrams of Fig. 3e and penguins with left coupling 
of gauge bosons to leptons: 



]TzM u ){(o*) 2 o^2Gx>M 

m 2 (Ol) 2 (B + B 1 )} 



M 



A 

£3,2 



2T(p2,Pz)S{p4,Pi) 



x D 



M 



£3 3 = 2T {P2,P3)S{p4,Pl) 



(C22) 



{u){{Of) 2 OfO z R 2C 00 {u) 
(of) 2 (Ol) 2 (B + Si)} (C23) 



x £ D v {u){(0t) 2 0f0 v R C(u) 

V=~/,Z° 

- (offiOlfiBv + B,)}. (C24) 
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Some important remarks: each diagram with a LFV 
and a LFC scalar line is described by the propagators of 
Eq. (9) and Eq. (10), so that the loop coefficients in the 
amplitudes are a sum of four integrals, while the ones 
with only the LFV line is a sum of two. The scalar two 
point function Bq and the tensor coefficients B\ , Coo that 
appear in the electroweak penguins are ultra-violet diver- 
gent, but the amplitudes are finite due the ortogonality 
of the slepton mixing matrix. 

Penguin diagrams with the exchange of the photon in 
the t or u channel are divergent for i, u — > 0. This di- 



vergence is cancelled by the graphs with external legs 
rcnormalization as required by gauge invariance. As ex- 
plained in Fig. 2, the ^-channel penguin diagrams where a 
scalar line is not dotted, contribute two times because the 
LFV propagator may appear once in both lines. The two 
amplitudes are equal because of the symmetry property 
of LoopTools form factors giving in this way a factor 
of two, that is necessary for the cancellation of the small 
t or u divergence. Finally, each amplitude gets a factor 
/f^yr = from the loop normalization convention. 
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